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Summary. We study the evolution of a small perturbation of the equilibrium of 
a totally asymmetric one-dimensional interacting system. The model we take as 
example is Hammersley's process as seen from a tagged particle, which can be 
viewed as a process of interacting positive- valued stick heights on the sites of Z. It 
is known that under Euler scaling (space and time scale n) the empirical stick profile 
obeys the Burgers equation. We refine this result in two ways: If the process starts 
close enough to equilibrium, then over times n v for 1 < v < 3, and up to errors 
that vanish in hydrodynamic scale, the dynamics merely translates the initial stick 
configuration. A time evolution for the perturbation is visible under a particular 
family of scalings: over times n u , 1 < v < 3/2, a perturbation of order n 1 ~ l/ from 
equilibrium follows the inviscid Burgers equation. The results for the stick model 
are derived from asymptotic results for tagged particles in Hammersley's process. 
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1. Introduction. A number of recent papers have sought various refinements to 
the basic hydrodynamic limits of interacting particle systems. One type of refine- 
ment is to add a small perturbation to the equilibrium of the process, and study the 
time evolution of the perturbation. For asymmetric exclusion in dimensions 3 and 
higher, Esposito, Marra and Yau (1994) proved that under diffusive scaling (time 
scale n 2 and space scale n) a perturbation of order n~ l follows a conservation law 
with a diffusion term. The backdrop of this result is the standard hydrodynamic 
limit of asymmetric processes, which leads to a conservation law without diffusion 
term under Euler scaling (time scale and space scale both n). A context for the 
Esposito et al. result is the search for microscopic interpretations of the Navier- 
Stokes equations. We refer the reader to p. 185-188 in Kipnis and Landim's (1999) 
monograph for a description of this program and further references. 

Our paper looks at the question of Esposito et al. in one dimension. We add a 
perturbation of order to the equilibrium, (3 > 0. The perturbation vanishes in 
the hydrodynamic limit n — > oo, and we study the effect of this perturbation under 
various time scales n v t, v > 1. We have two types of results: (1) For f3 G (0, 1/2), 
a hydrodynamic limit in the time scale n 1+/3 t shows that the perturbation obeys 
macroscopically the Burgers equation without diffusion term. (2) For u e [1, 3) and 
(3 close enough to 1, we show that the dynamics is simply a translation of the initial 
configuration, up to o(n) error terms. 

The most popular models for studies of the hydrodynamics of asymmetric sto- 
chastic dynamics are the exclusion and zero-range processes. Instead of these pro- 
cesses, we prove our results for the so-called stick process, which can also be regarded 
as Hammersley's process as seen from a tagged particle. This process has nonneg- 
ative variables r/(i) [stick heights] on the sites of Z that exchange pieces between 
each other. The stick process lacks some of the good properties of the exclusion or 
zero-range process: the state space is not compact, the rates are unbounded, and 
the amount of material that jumps is also unbounded. 

To make up for these complications the totally asymmetric one-dimensional pro- 
cess has a beautiful combinatorial structure uncovered by Aldous and Diaconis 

(1995) . This structure connects Hammersley's process and the stick model to the 
increasing sequences problem on planar Poisson points. A key ingredient of our 
proofs are sharp deviation estimates for the inceasing sequences problem from Kim 

(1996) , Seppaiainen (1998b), and Baik, Deift and Johansson (1999). 

We believe that the results of our paper hold also for totally asymmetric exclu- 
sion and zero-range processes. The basis for this conjecture is that these pro- 
cesses possess particle-level variational formulations that involve planar growth 
models, analogous to the increasing sequences connection of Hammersley's process 
[Seppaiainen (1998a,c)]. Johansson (1999) has shown that the limiting fluctuations 
for this growth model are the same as for the increasing sequences model. 

Organization of the paper. In section 2 we describe the stick model and state 
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the results mentioned above: Theorem 1 gives the translation, and Theorem 2 the 
hydrodynamic limit of the perturbation. Theorems 1 and 2 are corollaries of cor- 
responding Theorems 3 and 4 for tagged particles in Hammersley's process. These 
are stated in Section 3. The translation Theorem 3 for Hammersley's process is 
compared to a similar result of Ferrari and Fontes (1994) for asymmetric exclusion. 
Section 4 addresses briefly the rigorous construction of Hammersley's process and 
the stick process, and the connection with increasing sequences. Sections 5-7 con- 
tain the proofs: Section 5 contains lemmas, Section 6 the proof of Theorem 3, and 
Section 7 the proof of Theorem 4. 



2. The stick model and the results. Here is an informal description of the 
model. A rigorous construction will follow in Section 4. The state of the process is a 
configuration rj = (rj(i) : i G Z) where each rj(i) is a nonnegative real number. Think 
of r)(i) as the height of a vertical stick attached to site i G Z. At exponential rate 
equal to rj(i), the following event takes place: Pick a random quantity u uniformly 
distributed on [0, rj(i)]. Break off a piece of length u from the stick at i, and attach 
this piece to the stick at site i + 1. Thus if the neighboring stick lengths before the 
event were (r}(i),rj(i + 1)), then after the event they are (rj(i) — u, rj(i + 1) + u). 
These events happen at all sites i independently of each other. In the language of 
generators, this dynamics is expressed as 

(2-1) Lf( V ) = £ r lf(v u ' l ' t+1 ) - f(v)]du , 

where rj u,t,l+ = rj — u5i + u5i+i. 

In Seppalainen (1996) a Markov process rj(t) = (rj(i,t) : % G Z), t > 0, is 
constructed that operates according to this description. The state space of the 
process is 

(2.2) Y = \r} G [0, oo) z : lim A^" 2 V r}{i) = X 

k i=N J 

and the paths of the process are in the Skorohod space -D([0, oo), Y). Y is not 
closed in the product topology, but is given a stronger topology with a complete, 
separable metric. L of (2.1) is the generator of the process, in the sense that 

(2.3) E^[f( V (t))]-f( V ) = f E^[Lf( V (s))}ds 

Jo 



for all bounded continuous cylinder functions /on7 and all initial states rj G Y . E v 
stands for the expectation under the path measure of the process started at state 
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rj. Furthermore, the process has a 1-parameter family of invariant distributions, 
namely the i.i.d. exponential distributions on the variables (rj(i) : i G Z). 

We focus now on the hydrodynamic behavior of this process. The basic result 
[Seppaiainen (1996)] is that under Euler scaling the empirical stick profile obeys 
the Burgers equation. Suppose u(x,t), (x,t) G R x [0, oo), is an entropy solution 
of the Burgers equation 

(2.4) u t + (u 2 ) x = , u(x, 0) = u (x) 

with nonnegative initial data uq G L°°(R). Consider a sequence ?y n , n = 1, 2, 3, . . . , 
of stick processes, and assume that a law of large numbers is satisfied at time t = 0: 
for all a < b in R, 

\ nh \ rb 

(2.5) lim n _1 7 r] n (i,0) = / uo(x)dx in probability. 

i=[na]+l Ja 

The theorem is that the law of large numbers continues to hold at all later times 
t > 0: 



(2.6) 



lim n~ x } rj n (i,nt)= I u{x,t)dx in probability. 

™ , ^,1 J a 



i=[na] + l 



Euler scaling refers to the scaling in the above limit, where the ratio of macrosopic 
and microscopic units is the same n for both space and time. A macroscopic space 
interval (a, b] corresponds to approximately n(b — a) microscopic lattice sites, and 
macrosopic time t corresponds to microscopic time nt. The derivation of the hydro- 
dynamic limit (2.6) from the hypothesis (2.5) requires some technical assumptions, 
and the details can be found in Seppaiainen (1996). 

A trivial special case of the hydrodynamic limit is of course the case of a process in 
equilibrium: If the sticks are initially i.i.d. exponentially distributed with common 
expectation E[rj(i,0)] = q, then this situation persists, and the macroscopic profile 
is the constant u(x, t) = q. 

In the present paper we study the evolution of a small perturbation of the equi- 
librium. The initial macroscopic profile is 

(2.7) u {x) = q + n~ p v {x) 

where q > is the fixed equilibrium density, vq is a bounded measurable function 
on R, and f3 G (0, 1] is a parameter that we adjust to investigate different scalings. 
The function v is not assumed to take any particular sign, so to have nonnegative 
profiles we consider only n large enough to have q > n _ ^||wo||oo- For each n, the 
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initial stick configuration (?? n (i,0) : i G Z) is assumed to be in local equilibrium 
with macroscopic profile uq. Precisely speaking, our assumption is this: 

The variables (r] n (i,0) : i G Z) are mutually independent, 
exponentially distributed, and have expectations 

/■i/n 

E[r) n (i,0)] = q + n 1 '? v (x)dx . 

J(i-l)/n 

The perturbation of the expected density is taken to be 

n~P ■ {the average of vq ° ver the interval ((i — l)/n, i/n] }, 

instead of ■ {the point value vo(i/n)} because we are not making any regularity 
assumptions 011% A standing assumption is also that q > 0. In Section 4.1 we 
explain how the case q = is reduced to the basic hydrodynamic limit (2.5)-(2.6). 

Since \\uo — q\\oo — ► as n — > 00, the limit (2.6) is valid again with constant 
profile u(x, t) = q. To escape the regime of (2.6) we subtract the equilibrium density 
q and speed up time more, beyond the hydrodynamic scale nt. We introduce a 
second parameter v G [1, 00), and look at the evolution of the stick profile over 
times of order n u . The space scaling will be the same as in (2.5)-(2.6), so the 
lattice of sites scales as n~ 1 Z. Our object of study is the empirical stick profile 

52{v n (i,n"t)-q}6 i/n . 
iez 

In other words, we follow either integrals '^2 ie 2,{ r J n (h ri v t) —q}(f>(i/n) of compactly 
supported, continuous test functions 0, or equivalently the total stick mass in 
macroscopic intervals (x, y], Xll^Lxj+i V n (h nV i) ~ n l{y — x )- 

Let us derive an easy "benchmark" result against which we can compare later 
results. It is proved in Seppalainen (1996) that the stick process is attractive. This 
means that, if rj and £ are two initial states that satisfy 77 > ( [inequalities are 
interpreted coordinatewise, 77 (i) > ((i) for all % G Z], then it is possible to construct 
the processes 77 (t) and ((t) on a common probability space so that the inequality 
v(t) > C(t) holds at all times t > 0, almost surely. 

Fix ?i for the moment. Let C 1 and C 2 be stick processes in equilibrium, with 
expectations 

E[( 1 (i,i)]=q-n-P\\v \\ oo and E[( 2 {i,t)] = q + n~ 13 \\v \\ 00 ■ 

In other words, for each fixed time t and for r G {1, 2}, the stick heights (C r (i, t) : i G 
Z) are exponentially distributed i.i.d. random variables with expectations as above. 



(2.8) 
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At time t = 0, we can construct the initial configurations of all three processes C 1 , 
if 1 , and ( 2 on a single probability space so that 

0) < r] n (i, 0) < ( 2 (i, 0) for all i, a.s. 

To do this take an i.i.d. sequence of Exp(l) variables {Xi} and set 

C(i, 0) = E[( r (i, 0)]X { and r] n (i, 0) = E[n n (i, t)]Xi for r = 1, 2 and all i. 

We then construct all three processes on a common probability space so that 
CHM) < <(M) < C 2 (M) for all t and i, a.s. 

Let {-f^n} be an arbitrary sequence of integers, to be used as translations on the 
lattice. The construction gives these inequalities: 

[ny] 

{( 1 (K n + i,n"t) - E[{\i,n"t)]} - ([ny] - [na;] )n^||i;o||oo 

i=[nx] + l 
[ny] 

(2.9) < {v n (K n + i,n»t)-q} 

i=[nx] + l 
[ny] 

< {C 2 ( K n + h n v t) -E[( 2 (i, n't)] } + ([ny] - M)n^KHoo • 

i=[nx] + l 

A terminological convention: Throughout, we shall say 

(2.10) X = Y + o(n a ) a.s. 
as a shorthand for 

(2.11) lim n~ a \X-Y\ =0 a.s. 

The sums of the £-terms in (2.9) are almost surely ©(n 1 / 2 " 1 "* 5 ) for any 5 > because 
the variables are i.i.d. and have sufficient moments. 
Thus we get the result 

[ny] 

(2.12) V n (K n + i,n u t)=nq(y-x)+o(n [1/2Ml - f3]+s ) a.s. 

i=[nx] + l 

for any i/, (3 > 0, and arbitrarily small 5 > 0. The translation K n was included in 
anticipation of later results. Because we are speeding up time beyond the hydrody- 
namic scale, a certain translation will appear naturally. The goal of the remainder 
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of the paper is to improve on (2.12), by obtaining results that reveal how the per- 
turbation evolves in time, or have a smaller error term. 

As the last preparatory step, we construct the solution of the Burgers equation 
(2.4) by the Hopf-Lax formula. The perturbation vq(x) is now the initial data. 
Define Vq(x) by 

rv 

Vb(0) = and V (y) - V (x) = / v Q (z)dz for all x < y. 

J X 

V is a Lipschitz function with a bounded derivative a.e. For (x,t) G R x [0, oo), 
define V(x, 0) = V (x) and for t > 

(2.13) V(x,t)= inf !v (y) + hx-y) 2 

Then V is the unique viscosity solution of the Hamilton-Jacobi equation 

(2.14) Vt + ( Vx f = o, V(x,0) = V o (x). 

For each fixed t, V{- , t) is again a Lipschitz function, so it has derivative 
v = V x . This function v(x,t) is the unique entropy solution of (2.4) with initial 
data i>o- The reader can find a development of these p.d.e. results in Evans (1998). 

Now the results for the stick process. The most general result, valid for all 
scalings, does not identify any time evolution, only a translation of the initial sticks. 

Theorem 1. Assume that (3 > and v > 1. Let r] n (t) denote the stick process 
started from the initial configuration (2.8). Fix x < y in R and t > 0. Then, for 
any 5 > 0, the following asymptotic equality is valid almost surely as n — > oo: 

[ny] [ny] 

(2.15) V n ([2n iy qt]+i,n u t)= ^{h 0) + (n^- 2l3 ^^ +s ) . 

i=[nx] + l i=[nx] + l 

Why the translation \ln v qt\ appears naturally is explained in Section 4.2. The 
error exponent in the statement (2.15) satisfies 



(2.16) \v - 2(3] V [u/3] 



v - 2/3, v > 3/3 
i//3, z/<3/3 



If v < 3 and /3 > (z/ — 1)/2 the error in (2.15) is o{n) and so vanishes in the standard 
hydrodynamic scaling of (2.6). 

We find one family of scalings where the perturbation evolves according to the 
Burgers equation. For this to happen the perturbation has to be larger than n -1 / 2 . 
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Theorem 2. Suppose (3 G (0, 1/2) and set v = 1 + f3. Let r] n {t) denote the stick 
process started from the initial configuration (2.8). Let 4> be a compactly supported, 
continuous test function on R. Then almost surely 

(2.17) lim — l —, V{?7 n ([272 1+ V] + i,n 1+f3 t) -q}<j>{i/n) = [ (f)(x)v(x,t)dx . 

This result compares directly with Corollary 2.3 in Esposito et al. (1994), where 
the corresponding result is proved in dimensions d > 3 for an exclusion process. 
The deterministic limit (2.17) cannot be valid for (3 = 1/2 because in equilibrium 
this would be the central limit theorem scaling. 

Remark about construction. The almost sure results of Theorems 1 and 2 are proved 
for a special construction explained in Section 4. In this construction the processes 
rf 1 are defined on one common probability space, and the variables n n (i, t) are real- 
ized as interparticle distances of Hammersley's process. This is not the construction 
used in (2.9) that makes r\ attractive. Both theorems are proved by Borel-Cantelli 
arguments, and the probability estimates for the arguments are derived with the 
help of the special construction. But once derived, the estimates are valid in all 
constructions because they are statements about the distributions of the processes. 
Hence Theorems 1 and 2 are valid for any construction of the stick process. 

To compare (2.17) directly with (2.12), we can write it in the form 

[ny] 

il n {[2n 1+ Pqt] + i,n 1+ H) 

(2.18) i=[nx] + l 

ry 

= nq(y - x) + n 1_/? / v(z, t)dz + o(n 1_/3 ) . 

J X 

A comparison of the errors in (2.12), (2.15), and (2.18) reveals that for v > 1 + (3 
the easy result (2.12) in fact has the smallest error. For v = 1 + (3 (2.18) has 
the smallest error, while for v < 1 + (3 it can be either one of (2.12) and (2.15), 
depending on the exact relation between (3 and v. 

3. Asymptotics for a tagged particle in Hammersley's process. In Ham- 
mersley's process a countable collection of point particles evolves on R according 
to the following rule: if x < y are two locations of neighboring particles, then with 
rate equal to the distance y — x the particle at y jumps to a randomly (uniformly) 
chosen location in the interval (x, y). All particles execute jumps independently of 
each other. 

This evolution can be graphically constructed with a rate one homogeneous 
Poisson point process on the space-time plane R x (0, oo): Suppose (x, t) is a point 
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of the Poisson process. Then at time t the leftmost particle in [x, oo) jumps to x. If 
the leftmost particle already happened to be at x , or if there is no leftmost particle 
in [x, oo), no jump takes place. This latter case can happen if there are infinitely 
many particles in some bounded interval. 

There is an obvious connection between Hammersley's process and our stick pro- 
cess. We assume that we can label the particles by integers in an order-preserving 
way. Let z(i,t) denote the position of particle % at time t. The assumption is 

(3.1) z{% — l,t) < z(i,t) for all i and t. 

Suppose we have constructed the process z(t) = (z(i, t) : i e Z) that operates 
according to the description above. Define 

(3.2) rj(i,t) = z(i,t) - z(i - l,t) forzGZ. 

Then it is clear that r](t) evolves as our stick process. When particle z(i) jumps to 
the left, stick rj(i) donates a piece to stick rj(i + 1). In particle system jargon, the 
stick process is Hammersley's process as seen from a tagged particle. What this 
means is that knowing rj(t) and the evolution of one particle z(j, t) is equivalent to 
knowing the process z(t). The simultaneous construction of Hammersley's process 
and the stick process is discussed in Section 4. 

Assume that the initial sticks 7] n (0) = (i] n (i 7 0) : i G Z) that satisfy (2.8) have 
been defined on some probability space. Initial particle configurations z n (0) = 
(z n (i, 0) : i e Z) are defined on this same probability space by 

i 

^(0,0) = 0, z n (i,0) = ^?? n (j,0) forz>0, 

i=i 

o 

and z n (i,0)= ^?? n (j,0) for i < 0. 
j=i+i 

Thus z n (i, 0) is a sum of independent exponential random variables with uniformly 
bounded expectations, and 

(3.4) E[z n (i,0)] =qi + n 1 - p V {i/n). 

The processes {z n (t)} are then constructed together on one probability space 
where the initial configurations {z n (0)} and the space-time Poisson points are in- 
dependent. All processes z n (t) use the same realization of the Poisson points to 
construct the dynamics. This is not really necessary because our a.s. results come 
from Borel-Cantelli arguments. But in the proof it is convenient to work with a 
single Poisson process and the family {z n (0)} of initial configurations, instead of 
giving each process z n (t) its own space-time Poisson process. 



(3.3) 
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Theorem 3. Assume that (3 > and v > 1. Let z n (i,t) denote Hammersley's 
process started from the initial configuration described in (3.3) and (2.8). Fixx G R 
and t > 0. Then, for any 5 > 0, we have the following asymptotic equality almost 
surely as n — > oo; 

(3.5) 2 n (M + [2nV], Wt) = + z n ([na;], 0) + ( n W-WM»/3]+5j _ 

Theorem 1 is an immediate consequence of (3.2) and Theorem 3. 

Ferrari and Fontes (1994) proved a translation result of this type for the exclusion 
process. Suppose for the moment that the r/(z)'s are occupation variables of totally 
asymmetric 1-dimensional simple exclusion in equilibrium at density p. Then the 
jumps of the ^-variables correspond to the current of particles. Statement (1.5) in 
Theorem 1 of Ferrari and Fontes (1994) implies that, in the I? sense as n — > oo, 

(3.6) z(0, nt) = ntp 2 + z([nth(p)},0) + o(n 1/2 ) , 

where h(p) = 2p — 1. This can be compared with our result for Hammersley's 
process: with v = 1 and (3 > 1/3, (3.5) implies that 

(3.7) z n (0, nt) = ntq 2 + z n (-[2tqn], 0) + o(n 1/3+s ) . 

The error is smaller in (3.7) than in (3.6), but the Ferrari- Fontes result is valid for 
more general asymmetric exclusions, not only for totally asymmetric. 

Next a result with explicit time evolution. In one of the cases treated by the 
next theorem, we will assume that Vq(x) has asymptotic slopes in the sense that 
these limits exist: 

(3.8) vq(— oo) = lim Y^i^l anc i W+oo) = lim Y^i^l . 

When this is the case, we define the piecewise linear "asymptotic profile" 

{vo(— oo)x , x < 
0, x = 

v (+oo)x , x > 

and its evolution for t > by 

(3.10) V OQ (x,t)= inf {v oo (y,0) + hx-y) 2 
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Theorem 4. Assume v > 3/3, in addition to the basic assumption (3 > and v > 
1. Let z n {i,t) denote Hammersley's process started from the initial configuration 
z n (i,0) described in (3.3) and (2.8). Fix x G R and t > 0. TTien toe /iave £/ie 
following asymptotic equalities, each statement valid almost surely as n — > oo. 
Case i: z/ > 1 + /?. Assume that the limits in (3.8) exzst. T/ien 

(3.11) z n ([nx] + \2n v qt\,n v t) = n v tq 2 + nxq + n^V^O, t) + o(n"- 2f3 ) . 
Case 2: v = 1 + (3. Then 

(3.12) z n ([nx] + \2n v qt\,n v t) = n v tq 2 + nxq + n 1-/ V(:r, t) + o(n 1_/3 ) . 
Case 3: 1 < v < 1 + (3. Then for any 5 > 0, 

(3.13) z n ([nx] + [2n u qtln v t) = n v tq 2 + nxq + n^Voix) + o^^M^R^ _ 

Remarks. Recall again (2.10)-(2.11) for the precise meaning of the almost sure 
o(n a ) error terms. In Case 1, the term nxq in (3.11) may or may not be included 
in the error o{n u ~ 2 ^), depending on whether v > 1 + 2(3 or not. The statement 
(3.11) for Case 1 does not improve (2.12) because the error n v ~ 2 ^ is strictly larger 
than n 1_/3 in this case. Theorem 2 follows from Case 2 above. 

The remark about construction at end of Section 2 applies here too. The proofs 
of Theorems 3 and 4 are Borel-Cantelli arguments that depend on estimates of the 
distributions of the processes, and hence are valid in all constructions. 

The three cases reveal the effect of the time scale on the evolution of the pertur- 
bation: For fast times v > 1 + (3 we only see the asymptotic effect Voo(0, t) which 
is independent of the reference point x. For slow times v < 1 + (3 we only see the 
initial perturbation Vq(x). And exactly at v = we see the perturbation evolve 

according to the Burgers equation. 

It remains to prove Theorems 3 and 4. This proof uses on a special construction 
of Hammersley's process in terms of increasing sequences of the space-time Poisson 
points. 

4. Graphical construction and increasing sequences. Consider a planar, rate 
one, homogeneous Poisson point process. A sequence (xi,ti), (x2, £2), • • • , (x m , t m ) 
of Poisson points is increasing if 

x\ < x 2 < ■ ■ ■ < x m and t\ < t 2 < ■ ■ ■ < t m . 

For arbitrary (a,s), (b,t) on the plane, define the random variable L((a,s), (b,t)) 
as the maximal number of Poisson points on an increasing sequence contained in 
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the rectangle (a, b] x (s,t\. Abbreviate L(6, t) = L((0, 0), (6, t)) for the case where 
the lower left corner is the origin. 
An inverse to L is defined by 

(4.1) T((a, s),m, r) = inf [h > : L((a, s), (a + h, s + r)) > m} . 

In words: r((a, s), m, r) is the minimal horizontal distance h for which the rectangle 
(a, a + h] x (s, s + r] contains an increasing sequence of m points. Again abbreviate 
r(m,r) = r((0,0), m,r). 

These random variables satisfy laws of large numbers: 

1 la 2 

(4.2) lim -L(s6, st) = 2Vbt and lim ~r([sa], s£) = — a.s. 

S^OO S S^OO S 4:t 

The existence of the limits follows from the subadditive ergodic theorem. The exact 
values were first calculated by Vershik and Kerov (1977). 

In the previous section we suggested how to construct Hammersley's process 
with a rate one space-time Poisson point process. The rule was that Poisson point 
(x,t) pulls the leftmost particle in [x, oo) to the location x at time t. As usual in 
particle system contexts, constructing the process rigorously from this description, 
on the infinite real line, needs a proof. 

We can take an elegant way out with the help of the increasing paths. Assume 
given an initial configuration z(0) = (z(i, 0) : i G Z) that satisfies the ordering 
convention (3.1) for t = 0. Given a realization of the Poisson points, define 

(4.3) z(k, t) = inf \z(i, 0) + r((z(i, 0), 0), k - i, t) \ 

for all k G Z and t > 0. In words: The potential locations of z(k,t) are all points 
x such that the rectangle (z(i, 0), x] x (0, t] contains an increasing sequence of k — i 
Poisson points. Of these potential locations z(k, t) chooses the leftmost. 

If we permit — oo as a value for z(k,i), (4.3) defines a process z(t) = (z(k,t) : 
k G Z) that satisfies (3.1). To rule out the possibility of jumping to — oo in finite 
time, define the state space 

(4.4) Z = {z = (z(i)) G R z : z(i - 1) < z(i) for all i, and lim r 2 z(i) = 0} . 

One can check that if z(0) G Z, then almost surely the infimum in (4.3) is always 
attained at some finite i and z(t) G Z for all t. Homogeneity of the space-time 
Poisson point process then implies that (4.3) defines a time-homogeneous Markov 
process z{t) with state space Z. 
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Definitions (2.2) and (4.4) show that Y maps injectively into Z through equations 
(3.3), and Z back onto Y through equation (3.2). So given an initial stick config- 
uration (rj(i, 0) : i G Z) in Y, we define an initial particle configuration z(0) G Z 
as in (3.3), then define the process z(t) by (4.3), and finally use (3.2) to define the 
stick process r](t). 

This is convenient as a rigorous definition of the processes rj(t) and z(t), but it 
is not so obvious that the resulting dynamics follows our earlier descriptions. One 
can prove that when r](t) is defined this way, (2.3) is satisfied so the generator of 
r]{t) is L. All the facts mentioned here can be found in Sections 3-5 in Seppalainen 
(1996). 

We can also argue from definition (4.3) that if (x,t) is a space-time Poisson 
point, then at time t the leftmost particle in [x, oo) is at x, if such a particle exists. 
Suppose not, so that for some k, z(k — 1, t) < x < z{k, t). Pick i < k — 1 so that 

z(k -l,t) = z(i, 0) + r((z(i, 0), 0), k - 1 - i,t) . 

Barring the null event that space-time Poisson points can lie on the same horizontal 
line, there must be an increasing sequence of k — 1 — i Poisson points from (z(i, 0), 0) 
to a point (y, s) such that y = z(k — l,t) < x and s < t. [In the extreme case 
i = k — 1, this sequence is empty, and s = 0.] Consequently we can append the new 
point (x,t) to this increasing sequence to produce a sequence of k — i points from 
(z(i,0),0) to (x,t). Then 

z(k, t) < z(i, 0) + r((z(i, 0), 0), k - i, t) < x , 

contradicting x < z{k,t). 

The remainder of the paper proves Theorems 3 and 4 through definition (4.3). 
The construction of the family of processes {z n (t)} is the following. There is a single 
probability space (0,JF, P) on which are defined the initial locations {z n (0)} and, 
independently of them, the space-time Poisson point process. On this probability 
space define the random variables 

(4.5) r n (i,m,t) = T((z n (i,0),0),m,t) . 
Then, following (4.3), the processes {z n (t)} are defined by 

(4.6) z n (k,t)= inf \z n (i,0) + r n (i,k-i,t) 1. 

i:i<k I J 

Our arguments use distributional bounds on the initial locations z n (i,0) and the 
variables _T n (i, m, t) . In distribution r n (z, m, t) is equal to T(m, t), so we can ignore 
the indices n and i and switch to T(m, t) as soon as only distributional properties 
are studied. 

We close this section with two comments about matters that came up in Section 

2. 
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4-1 The case q = 0. Our results in Sections 2 and 3 are for the case where the 
fixed equilibrium density q is stricly positive. Here we show how the case q = 
reduces to the standard hydrodynamic setting, through a space-time scaling of the 
graphical picture. 

Suppose q = 0, and let the initial configurations i] n (0) and z n (0) be as in (2.8) 
and (3.3). Define another initial particle configuration by z n (0) = n l3 z n (0). Con- 
struct the two processes z n (t) and z n (t) by formula (4.3), with these Poisson pro- 
cesses: for z n (t) take a realization IT of the rate one space-time Poisson points, 
and for z n (t) use the space-time points IT obtained by mapping the points of IT 
by (x,t) I— > (n^x,n _/3 t). By the scaling properties of Poisson processes, IT is again 
a rate one homogeneous Poisson point process. But notice that the difference be- 
tween the evolutions z n (t) and z n (t) is merely a stretching and shrinking of the 
space and time axes: z n (i,t) = n~^z n (i,n~^t). By the standard hydrodynamic 
result n~ 1 z n ([nx] 7 nt) — > V(x,t), so we get that n^~ x z n {]nx\ 1 n 1+ @t) — > V(x,t). 
This is Case 2 of Theorem 4. And similarly, Theorem 2 holds for the stick model. 

4-2 The translation \2n v qt\. We advance here some explanation for the spatial 
translation [2n v qt] in our theorems. 

First consider the variational formula (4.3). Suppose that the process is in equi- 
librium with E[r](i,t)} = q. (4.2)-(4.3) give 

(4.7) z(0, n u t) = inf iqy + + [fluctuations] 1 . 

y<o i An v t J 

Neglecting fluctuations, the infimum is attained at y = —2n v qt. So roughly speak- 
ing 

^0,n^) = ^(-[2nV],0) + r((^(-[2nV],0),0),[2nV],^) 
+ [fluctuations] . 

As a sum of independents the term z(— [2n u qt], 0) has fluctuations of order n^/ 2 , 
while the T-term has fluctuations of order ?W 3 (Lemma 5.2). It is advantageous 
to move the translation [2n u qt] to the left-hand side of (4.8), so that we study 
the dynamics of z($.n u qt\,n u t). Then the minimizer in (4.7) is y = 0, and we get 
smaller fluctuations on the right-hand side of (4.8). 

Alternatively, we can look at the macroscopic equation to find the right scaling 
and translation for nontrivial dynamics. Suppose first that u(x, t) = q + n~ l3 p(x, t) 
satisfies the Burgers equation (2.4). Then p(x,t) satisfies the equation 

p t + 2qp x + n- f3 (p 2 ) x = 0. 

In the limit n — > oo this gives p t + 2qp x = which is solved by a spatial translation 
p(x, t) = po(x — 2qt). To get nontrivial dynamics we speed up time and set w(x, t) = 
p(x,n l3 t) that satisfies 

w t + 2n l3 qw x + (w 2 ) x = 0. 
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To eliminate the n^-term let v(x,t) = w(x + 2n f3 qt,t). The function v then 
solves Burgers equation vt + (v 2 ) x = again. Working backwards, v(x,t) = 
p{x + 2n p qt,n p t). 

To see how v(x,t) should arise microscopically, start with the "hydrodynamic 
heuristic" u(x,t) « (2ne)~ 1 Yl,\i\<ne v([ nx ] + h n t)- From this, 



POM) ~ Y £ ' ^i^F ^2 {v([nx]+i,nt) - q}, 

\i\<ne 



and then 



v(x, f ) « 1 . -1^ £ [v([nx] + [2n 1+ V] + i, n^H) - q] 

\i\<ne 



This is exactly what Theorem 2 states, with the right translation [2n 1+ ^qt] again. 



5. Auxiliary lemmas. Throughout the proofs, we use symbols C, C\, C2, etc, 
for constants independent of the important indices of the proof (such as m, n, i, or 
j). The values of C, Ci, C 2 , ... may change freely from one inequality to the next. 

We start with an inequality for bounding the initial locations z n (i, 0). 



Lemma 5.1. Suppose {Xi} are independent exponentially distributed random vari- 
ables with expectations E[Xi] = qi e [0, b] where b is a finite constant. Then for all 
e G (0,1/2) there is a finite constant C = C(b,e) > such that for large enough 
m E N, 



(5.1) 



s m m ^ 

Pi J2 X *~J2^ >ern 1 / 2+£ \<2exp(-Cm 2e ) 
< i=i i=i > 
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Proof. The standard exponential Chebyshev argument. Let t G (0, 1/6). 

, m m 

H 5>i > 5> + em 1/2+e 

^i=l i=l 

s m n in 

< expi -tJ2<li- tem^ 2+e \ J] E [e tx *] 

^ i=\ ' i=l 

s m m \ 

= expi -tj^qt- tem^ 2+e - ^ log(l - tq % ) \ 
^ i=i i=i ' 

= expj-tem 1 /^ + g (?*?- + 0(t 3 )) } 

< exp{-fem 1/2+e + mb 2 t 2 /2 + 0(t 3 m)} 

[ choose t = 6~ 2 £m e_1/2 ] 
= exp{-£ 2 &- 2 m 2 72 + 0(m 3e " 1/2 )} 

< exp(-Cm 2e ) . 

The last step is valid for e < 1/2. In the second equality we expanded log(l — 
tqi) = —tqi — q 2 t 2 /2 + 0(t 3 ) where the O-term is uniform over % because of the 
uniform bound qi < b. The expansion is valid because t = b~ 2 em e ~ 1 / 2 can be 
made arbitrarily small by restricting m to be large. The same argument with some 
sign changes proves also the other inequality. □ 

Next bounds on the fluctuations of the increasing paths. 

Lemma 5.2. Suppose a, s and h are positive real numbers. 

(a) For x > 2, define 

(5.2) I(x) = 2xcosh _1 (x/2) -2y/x 2 -4 . 
Then for a < hs < a 2 / 4, 

(5.3) P{r([a], s) < ^ - h} < expj-^vV-4/is j( 2 + ^) } . 

When x = hs/a 2 is small, we can use the expansion 

(5.4) 1(2 + x) > Cx 3/2 . 

(b) There are fixed positive constants B , B\, do, Co and C\ such that if a > B 
and Sia 4 / 3 < hs < doa 2 , then 

( a 2 } ( s 3 h 3 

(5.5) P\T{[a],8)>- + h\<C expj -d 
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Proof. Part (a). The random variables L(s, s) are superadditive is the sense that, 
for any < s < t, 

L((0, 0), (s, s)) + L((s, s), (£, £)) < L((0, 0), (£, *)) a.s. 

It follows that there exists a function 7(x) such that 

(5.6) sup - log P{L(s, s) > sx} = lim - log P{L(s, s) > sx} = I(x) . 

s>0 S s^oo s 

Since s _1 L(s, s) — > 2 as s — > oo, /(x) = for a; < 2. Kim (1996) proved that I(x) 
is bounded below by the expression in (5.2), and Seppalainen (1998b) that this 

expression equals I(x). By (5.6) and the observation that L(a, b) = L(v / a6", v^aF), 

p{r([a], a ) < ^ - 4 = P { L (£ - s ) - [a] 

2[al 



< exp < — v a 2 — 4/is I 



Va 2 - 4/is 
The argument of /(•) is estimated below by 

2[«] . 2(1-1^ ^ 1^ ( faU., fa 



Va 2 - 4/is Vl - 4/isa" 2 \ a J \ a? J a 
provided hs > a. 

Part (b) is a consequence of case 4 of Lemma 7.1 in Baik et al. (1999). We check 
the assumptions of that Lemma. First we express the probability (5.5) in terms of 
L, then convert it to the n ( Annotation of Baik et al. (1999): 

p{r([fl], 8 ) > ^ + h} = p{l(£ + < [a] } 

According to case 4 of Lemma 7.1 in Baik et al., we can bound 
(5.7) [aM ^ + fa) < C exp(C 1 t 3 ) 

with t defined by the equation 
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provided 



(5.9) 



1 + 



M 7 



< 



Va 2 + Ahs 



< 1 + S 6 , 



2 1 /3[ a ]2/3 



[a] 



where M7 and ^ are certain positive constants that appear in the development of 
Baik et al. 

The first inequality in (5.9) is equivalent to 



Thus the first inequality in (5.9) is satisfied if £?ia 4 / 3 < hs < 2a 2 for a large enough 
constant B\. 

For the second inequality in (5.9) observe that 



which is < 1 + $6 provided a is large enough and hs < cfoa 2 for a small enough do. 

We have verified the conditions of case 4 of Lemma 7.1 in Baik et al. Again 
because (5.11) is below the right-hand side of (5.10), we see that t defined by (5.8) 
satisfies —t > 2 1 / 3 /isa -4 / 3 , so inequality (5.7) becomes (5.5). □ 



Remark 5.1. In our typical application of Lemma 5.2, a and s are of the same large 
order m, and h is of the order m 1 / 3+e . Then the bound in (5.3) is C\ exp(— C2m 3£ / 2 ) 
and in (5.5) C\ exp(— C2m 3e ). 

Remark 5.2. For the growth model associated with the exclusion and zero-range 
processes, the lower tail estimate (5.3) is available in Seppalainen (1998c) and in 
Johansson (1999). But the upper tail estimate (5.5) has not been derived at the 
time of writing this paper. 



(5.10) 




Provided hs < 2a 2 , the right-hand side above is bounded below by 



(5.11) 
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Lemma 5.3. For e > 0, a > 0, (3 > 0, and v > 1, define a deterministic quantity 
R n by 

(5.12) R n = n v tq 2 + nxq + n^KIUM + a(n" /3+£ + n 1/2+£ ) . 
Then there are finite constants Ci > such that 

(5.13) P{z n ([nx] + [2n v qt],n v t) > R n ) < d exp(-C 2 n 2e ) 
for all n. 

Proof. By the variational formula (4.6) and Lemmas 5.1 and 5.2, 

P{z n ([nx] + [2n u qt],n v t) > R n ) 

< P{z n ([nx],0) > nxq + U 1 - 13 \\vqWoo\x\ + an 1/2+e ] 

+ P{r n ([nx], [2n v qt],n v t) > n v tq 2 + an v ^ +e ) 

< d exp(-C 2 n 2£ ) + Ci exp(-C 2 ?i 3e ) . □ 

The main lemma of this section reduces the range of indices that need to be 
considered in the variational formula (4.6). 

Lemma 5.4. Suppose (3 > and v > 1, and let £ be any number that satisfies 
£ > v — (3 and £ > 2v/3. Then if b > is large enough, the following holds with 
probability one: for large enough n, 

(5.14) z n ([nx] + [2n v qt],n v t) 

= min{z n (z,0) + T n (z, [nx] + [2n v qt] -i,n v t) : \i - [nx] | < [6n ? ] }. 

Furthermore, 

(5.15) P i ( 5 - 14 ) fails forn } < oo. 

n 

Remark 5.3. The distinction between the weak inequality £ > v — (3 and the strict 
inequality £ > 2i//3 in the hypothesis is significant. It arises in the proof of this 
Lemma 5.4 and influences the error terms of our theorems. 

Proof of Lemma 5.4- We shall prove (5.14). The reader can accumulate the esti- 
mates as we proceed and observe that (5.15) is also true. 
By Lemma 5.3 and Borel-Cantelli, we may suppose that 



z n ([nx] + [2n v qt],n v t) < R 
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at the expense of discarding an event of probability zero and by taking n large 
enough. Fix 5q G (0, 1) and take e small enough in the definition (5.12) of R n 
so that z//3 + e < v and 1/2 + e < 1. Then, since \z n {[{\ + 5 )n"tq] + [nx],0)\ 
is a sum of |[(1 + 5o)n v tq] + [nx]\ independent exponential random variables with 
expectations in the range q ± n _ ^||i>o||oo, basic large deviation estimates show that 



(5.16) ^P{z n ([(l + 5 )n v tq] + [nx],0) < R n ) < oo . 

n 

Consequently we may also assume 

(5.17) z n ([(l + 6 )n»tq] + [nx},0) > R n . 

It then follows from the variational formula (4.6) that, almost surely, 

(5.18) z n ([nx\ + [2n u qi\,n v t) 

= M{z n (j, 0) + r n (J, H + [2n V] " h n v t) : 

j < [(1 + 6 )n v tq] + [nx]} 

for large enough n. 

To conclude the proof we shall show that indices j outside the range \ j — [nx] \ < 
[bn^] cannot give the infimum in (5.18). This will follow from showing that, almost 
surely, 

z n ([nx] + [6n € ] +i,0) + r n ([nx] + [bn*] +i, [2n v qt\ - [bn^] -%,n v t) 
' - 1 >z n ([nx],0) + r n ([nx],[2n v qt],n v t) 

for all < i < [(1 + 6 )n v tq], and that 

z n ([nx] - [but] - i, 0) + r n ([nx] - [bn^] - i, [2n v qt] + [bn^] + i, n v t) 
{d ' 2U) > z n ([nx],0) +r n ([nx], [2n v qt],n v t) 

for all % > 0. The upper bound j < [(1 + 5$)n u tq] + [nx] in (5.18) permitted us 
to restrict the range of % to % < [(1 + ^o)^.^^^] in (5.19). The benefit is that the 
argument [2n u qt] — [bn^] — i of _T n in (5.19) is of order n u throughout the range of z, 
which makes the estimation easier because there is no need for separate arguments 
for values of smaller order. 

We estimate various terms separately in three sublemmas. 
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Sublemma 5.1. For any fixed b > and 5 > 0, the following statements hold 
almost surely: for all large enough n, 

(5.21) r n ([nx] + [bu^] + i, [2n u qt] - [bn^] - i, n u t) 

> (2n"qt -bn^-i) 2 - 8n" /3+S " 
4n»t K ' 

for allO<i< [(1 + 6 )n v tq], and 

(5.22) r n ([nx] - [but] - i, [2n v qt] + [bn^] + i, n v t) 

> J— (2nV + bnt + i) 2 - 5(n»/ 3+s » + 

An v t v ' v ' 

for all i > 0. 

Proof of Sublemma 5.1. We shall prove (5.22) and leave (5.21) to the reader. Their 
difference is that in (5.22), due to the unbounded range of i, we need % explicitly in 
the estimates and sum over % > in the end. (5.21) is easier because one estimate 
uniformly over i is sufficient. 

Let A n denote the event that (5.22) holds for all i > 0. Our goal is to prove 

(5.23) £P(^)<oo 

n 

so that by Borel-Cantelli A n happens for all large enough n, a.s. For fixed n and i, 
the event that (5.22) fails has the same probability as the event 

(5.24) T([2n v qt] + [bn 6 ]+i,n v t) < t^— ( 2 ™V + bn? + if - 5{n v/3+5v + i 1/3+s ) . 

By shrinking b slightly we can discard the integer parts. We bound the probability 
of the event (5.24) by (5.3). Now 

1/2 



\/ a 2 — 4hs 



(2nV + ^ + i) - 4^(n 4l/ / 3+5l/ + z^+V) 



>Ci(n" + z) 

and 

/, „ rt 'iv/3+6v i Al/3+5 rt i/ 

gt n v + i 

uniformly over % > 0, for all large enough n. Thus the probability of the event 
(5.24) is at most 

exp[-C(n v + i)A2 + C in - 2 ^ s+s n] < exp(-Cn 3 ^/ 2 - Cm*/ 2 ) , 
where we applied the expansion (5.4). Summing this over z > 0, we get that 

P(A c n ) < ^exp(-Cn 3 ^/ 2 - Cin-» +s5 »/ 2 ) < Cm^ 1 " 35 / 2 ) exp(-Cn 3 ^/ 2 ) . 

i>0 

This last expression is summable over n, so (5.23) happens. □ 



22 



Sublemma 5.2. For any fixed 5 > 0, this holds almost surely: for all large enough 
n, 



Proof of Sublemma 5.2. Deviation bound (5.5) and Borel-Cantelli. □ 

Sublemma 5.3. For any fixed b > and 5 > ; the following statements hold 
almost surely: for all large enough n, 



for allO<i< [(1 + 6 )n"tq], and 

(5.27) z n ([nx], 0) - z n ([nx] - [bn^\ - i, 0) 

< ( bn Z + i)( q + n-"||t;o||oo) +S(n^ 2 + 5 + i l ' 2+s ) 



for all i > 0. 

Proof of Sublemma 5.3. This lemma is a consequence of Borel-Cantelli and the 
distribution of the increments z n (j + l, 0)—z n (j, 0). We prove (5.27). The argument 
for (5.26) is the same. 

We can realize the initial locations z n (j, 0) so that the inequalities 



(5.28) (q-n-P\\v \\ oo )X j <z n (j,0)-z n (j-l,0) < (g + n-*> ||oo)*j 



are valid, where {Xj} are i.i.d. exponential random variables with expectation 
E[Xj] = 1. For fixed n and i, the probability that (5.27) fails is bounded above by 



(5.25) 



r n ([nx], [2n"qt],n v t) < n v tq 2 + 6n 1//3+5l/ . 



(5.26) z n ([nx], 0) - z n ([nx] + [bn*] + i, 0) 



< -{bri* + i)(q- rr"|M|oo) +d{n^ 2+s + t 1 ' 2 ^) 



(5.29) 




This probability is bounded above by 



(5.30) 
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where k(x) = x — 1 — logx is the Cramer rate function for the Exp(l)-distribution. 
In case the reader is used to thinking of large deviation rate functions only asymp- 
totically valid, note that inequality P{^2™ Xj > ma) < exp{— mn(a)} for a > 1, 
and its lower tail counterpart, are valid for finite m due to the supermultiplicativity 
P{Y!i m X 3 > (/ + m)a)) > P(Ei Xj > la)) ■ P(£? *j > ma)). 

For small x we have the quadratic lower bound k(1 + x) > Cx 2 , so (5.30) is 
further bounded above by 

f r)?+ 2<5 4- ■i 1 + 2S 

(5.31) exp \-C 



bn^ + i 

Summing the quantities in (5.31) over i > 0, we bound the probability that, for a 
fixed n, (5.27) fails for some i > 0, by 

£+25 , ,-1+2(5 



i>0 k J 



< exp(-Cn 2<5 ) + exp(-Cz 2<5 ) . 

0<i<n« i>n« 

The last line above is summable over n. Hence by Borel-Cantelli, it is almost surely 
true that for large enough n, (5.27) holds for all i > 0. □ 

We return to complete the proof of Lemma 5.4. By (5.21), (5.25), and (5.26), 
inequality (5.19) will hold for large enough n if we can show that 

(5.32) ^ (2n V -bn^-if - n u tq 2 - 25n u/3+5v 

> + i)(q- n-^|KHoo) + S(nt /2+S + t 1/2+S ) 
holds for all < i < [(1 + 6 )n"tq]. (5.32) simplifies to 

(5.33) h l n 2 ^ + -in^ v + — - 25n»/ 3+s » 
v ' At 2t An v t 

> fell^Hoon^ + \\voWooin-P + 5(n^ 2+s + z 1 / 2 ^) . 

Now suppose £ > v — (3 so that n 2 ^~ v > n^~@ and n^~ v > . Then (5.33) follows 
from 



>S(2n^ 3+5 " + n^ 2+s + i 1 / 2+s ) 
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Now choose 6 > 4t||v ||oo so that the two coefficients in parentheses on the left-hand 
side are positive and large. The condition £ > 2v/3 is exactly what is needed to 
have n 2 ^~ u > n u ^ s+Sl/ + n^ 2+s for all large enough n, if 5 > is small enough. The 
i-term on the right-hand side is controlled by the observation 

5i l ' 2+s <C(n 2 ^ +in*- v ) 

for all % > 0, provided n is large enough. 

The argument for (5.20) goes exactly the same way. This completes the proof of 
Lemma 5.4. □ 



6. Proof of Theorem 3. Let £ satisfy £ > 2v/3 and v — (3 < £, < v — £ for some 
small e > 0. Let M < oo be a large finite constant, to be chosen later. Define the 
events 

A n = { \z n { [nx] + [2n V] , n v t) - z n ( [nx] , 0) - n u tq 2 \ > 2Mn 2i ~ v } . 
By Borel-Cantelli, Theorem 3 will follow from proving the summability 

(6.1) ^P(A n )<oo. 

n 

To see this, compare 2£ — v with the exponent in the error of (3.5): In the case 
v < 3/3, take £ = 2z//3 + 5/3, so that 2£ — v < u/3 + 5. In the case v > 3(3, we can 
take £ = v - (3 so that 2£ - v = v - 2(3 < v - 2(3 + 5. 

Now to prove (6.1). By Lemma 5.4, for large enough n it is the case that 

z n ( [nx] + [2n u qt] , n v t) - z n ( [nx] , 0) - n u tq 2 

(6.2) = , ,inf U n ([nx] - i, 0) - z n ([nx],0) + qi 

+ r n ([nx] - i, [2n u qt] + i, n v t) - n v tq 2 - qi) . 
The probability P(A n ) is bounded above by 

p { \z n ([nx] - i, 0) - z n ([nx], 0) + qi\ > Mn 2i ~ v ) 

\i\<bnt 

(6.3) + P{ \T([2n v qt] + i, n v t) - n v tq 2 - qi\ > Mn 2i ~ v } 

\i\<bni 

+ P{ (6.2) fails for n } . 
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To bound the first probability in (6.3), apply (5.28) to get, for each i, 
P{ \z n ([nx] -i,0)- z n ([nx], 0) + qi\ > Mn 2i ~ v ) 

r w i 

< P\ J](g - n^WM^Xj < q\i\ - Mn 2 ^ \ 

r M i 

+ P\ + n ~ fi \\ v o\\oo)Xj > q\i\ + Mn^~ v \ 

r w i 

< J2 X J ^ \i\ + C 1 n- p \i\-M 1 n 2 ^- v V 

r w i 

+ P<^X j > |z| -Cin _/3 |i| + Mm 2 ^ L 

In the second inequality, C\ is a new constant, and M\ = M/(q + 1) accounts for 
the effect of dividing M by (g ± « n|t>o||oo) when n is large enough. 

Since < 6n^ _/3 and £ > z/ — f3 =>- 2£ — v > £ — (3, by choosing M 

large enough at the outset we have a further constant M 2 > such that M\n 1 ^~ v 
— C\n~^\i\ > M 2 n 2 ^~ v . Next, apply the large deviation rate function k for Exp(l)- 
variables as in (5.29)-(5.30). The new upper bound becomes 

P^Xj < \i\ -M 2 n^- v \ +p!.J2Xj > \i\+M 2 n 2 t-A 

< exp{ — |z|«(l - M 2 n 2i ~ v \i\~ 1 )} 

(6.4) + exp{-|i|«(l + M 2 n 2 ^~ u \i\~ 1 )} 

< exp jiWsn 2 ^ + |i| log(l - M 2 n 2 ^|z| _1 ) j 

+ exp j-M 2 n 2 ^ + |i| log(l + A^ra 2 * - "!*! -1 ) j . 

Check that the functions (1/x) log(l ± x) are maximized by taking x > as small 
as possible. Thus we get an upper bound for (6.4) by replacing \i\ by its upper 
bound bn^. Expanding the log then gives the upper bound 2 exp(— Cn^~ 2v ) for 

(6.4) [here the assumption £ < v becomes useful]. Tracing backwards, we conclude 
that 

(6.5) [ the first sum in (6.3) ] < Cm c exp(-Cn 3C " 2iy ) . 
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(6.6) 



Now to bound the second probability in (6.3). Again because \i\ = O(n^), 
P{ \T([2n"qt] + i, n v t) - n v tq 2 - qi\ > Mn 2i ~ v } 

| T([2n v qt] +i,n u t) - -^—(2n v qt + if > Mm^j , 



< P< 



for a constant Mi > 0, provided M was chosen large enough. Since |i| < Cn^ < 
Cn u ~ e , Lemma 5.2 implies that the probability in (6.6) is at most C\ exp(— C2n 3 ^~ 2u ) 
It follows that the bound in (6.5) works also for the second sum in (6.3). Thus the 
sum in (6.1) has the following bound: 

^P{A n ) < exp(-Cn^- 2v ) + ^P{(6.2) fails for n) < oo . 



The summability is a consequence of the assumption £ > 2^/3 and (5.15). (6.1) 
holds, and we have proved Theorem 3. 



7. Proof of Theorem 4. The proofs for the different cases are Borel-Cantelli 
estimates for the distribution of the random variable 2 n ([nx] + [2n !/ gt], n v t). For the 
sake of readability, we do not formulate explicit probability estimates and instead 
write statements of the type (2.10)-(2.11). Behind each a.s. error estimate is a 
summable deviation probability, as the reader can verify from the arguments. 

Case 3 can be proved quickly from Theorem 3: in (3.5) replace the term z n ([nx], 0) 
by its expectation (3.4) plus fluctuation o{n l / 2+s ). We concentrate on proving Cases 
1 and 2. 

Proof of Theorem 4, Case 1. Assuming v > 1 + /?, the goal is to show that 
(7.1) lim n w - v \z n {[nx] + [2n"qt],n v t) -n v tq 2 -nxq] = V^O,*). 

n^oo I J 

Let y be a number that achieves the infimum in (3.10) for x = 0. Set i = \n v ~^y\ 
in the expression inside the braces on the right-hand side of (5.14). For large n, we 
get the upper bound 

z n ([nx] + [2n u qt],n u t) 

< z n ([n^y},0) + r n ([n^y], [nx] + [2n»qt] - [n^y],n"t) 

< n»- p yq + n 1 "^™" V^D 



(7.2) 



v 2 

+ ™ u tq 2 + ^-n"~ 213 + nxq - n u ~ p yq + o{n v ~ 2(i ) 
< n v tq 2 + nxq + n u ~ 2p ^-^V^'^y) + + o 



(n 
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The following steps were taken above: For small e > 0, it is almost surely true that, 
for large enough n, 

z n ([n»-Py},0) 

< E{z n {[n u - l3 y\^)}+en^-^' 2+£ 

(7 ' 3) = [n v -< 3 y)q + n 1 "^™" V _/3 ?/]) + en^~^l 2+e 

< n'-Pyq + n 1 - f3 V (n u - f3 - 1 y) + o{n v - 2<3 ) . 

The first step above is by Lemma 5.1. v > 3/3 guarantees that n^ _/3 ^ 2+e = 
o(n iy ~ 2f3 ) if e > is small enough. 

Similarly by (5.5) for large enough n, 

r n ([n v - p y], [nx] + [2n u qt] - [n v - p y],n v t) 

(7.4) < + [ 2n V] " [n^Py]) 2 +en u/3+£ 

v 2 

< n v tq 2 + y —n v ~ 2(i + nxq - n v ~ p yq + d + C 2 n 1 ~ f3 + en u / 3+£ . 

The term C\ + C 2 n 1 ~ (3 accounts for terms left out after expanding the square and 
for removal of integer parts [•]. The assumptions v > 1 + (3 and v > 3/3 guarantee 
that Ci + C 2 n 1 ~ f3 + e?W 3+e = o{n v ~ 2f3 ) if e > is small enough. 

(7.3) and (7.4) justify the validity of (7.2) for large enough n, almost surely. Now 
we can prove one half of (7.1): 



(7.5) 



limsupn 2/3 - !/ L n ([nx] + [2n"qt], n v t) - n v tq 2 - nxq \ 

n — »oo l_ J 

< lim^supjn 1 -^^-^) + = VoofaO) + ^ = V^t) . 



The second last equality follows from 

(7.6) lim m- 1 V (my) = V oo (y,0). 

m^oo 

It remains to bound the liminf in (7.1) from below. By the assumptions v > 3/3 
and v > 1 + /3, we can choose a number g that satisfies 

(7.7) 1 - (3 < g < v - 2/3 , g > u/3 > (3 , and q > (v - (3)/2 . 
Define a sequence of deterministic numbers by 

r n = n v tq 2 + nxq + min < n 1 ~ f3 Vo(i/n) H > — 2n e . 

\i\<bn"-e\ An v t) 



2<S 



Lemma 7.1. Almost surely, the inequality z n ([nx] + [2n u qt],n"t) > r n holds for 
large enough ra. 

Before proving the lemma, let us use it to finish the proof of Case 1 of Theorem 

4. 

liminf n 2p ~ u \z n ([nx] + [2n v qt],n v t) - n v tq 2 - nxql 

> liminf min (n 1 ~ t/+/3 V r (i/n) + ; * ofl 1 

[ change of variable i = n u ~ f3 y ] 

> liminf inf ( n 1 - !y+/ VoK~ 1 ~ /3 y) + —\ 
~ n^oo \ y \< b \ uv yj At) 

> ^oo(0,t). 

To check the last inequality, let rij be a subsequence along which the liminf n ^ oc is 
realized. For each j pick y n . that realizes the infimum, pass to a further convergent 
subsequence y Uj — > y, and now consider different cases: If y Uj stays bounded away 
from zero, it follows from (7.6) that 



(0,t) 



(7.9) lunj^-^Vb^- 1 -^.) + ^} = V oo (y,0) + y 2 /4t > V a 
And if y = 0, Lipschitz continuity of Vb gives 

(^"^^(nj -1- ^^-)! < ||v ||oo|2/n 3 | — ► 0, 

so in this case too the limit in (7.9) is 14o(?/7 0) + y 2 /4t = 0. 

(7.5) and (7.8) together imply (7.1), and thereby prove Case 1 of Theorem. 
Before moving on to Case 2, we check Lemma 7.1: 

Proof of Lemma 7. 1 . Abbreviate temporarily 

(7.10) Z n = mm{z n (i, 0) + T n (z, [rax] + [2n v qt] - i, n v t) : |z| < bn v ~ p }. 

The assumption v > 3/3 permits us to set £ = v — f3 in Lemma 5.4, so z n ([ni] + 
[2n"qt], n u t) — Z n for large enough ra. The difference between |i| < hn v ~^ in (7.10) 
and |i — [nx]\ < hn v ~^ in (5.14) is irrelevant now because v — (3 > 1 and we can 
always increase h. To prove Lemma 7.1, we show that 

(7.11) Z n > r n holds for large enough ra. 
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The complementary probability P{Z n < r n } is bounded above by the sum 

P{z n (i,0) < qi + r^-^V^i/n) - n Q ) 



i\<bn 



v-(3 



(7.12) + E P { rn (^ M + [Wqt] - i, n v t) 



i 2 

< n v tq 2 + nxq — qi + n e 



In the first sum above the term for i = vanishes because by construction z n (0, 0) 
= Vo(0) = with probability one. We bound the sum over 1 < i < bn v ~^ and leave 
the matching argument for negative i's to the reader. Let {X,} be as in (5.28). 
First split the sum. 

P{z n (i,0) < qi + r^-PVoii/n) - n e } 

\<i<bn v -f 3 

(7.13) < P {Z)^- w "^INI|oo)^ < qi + n^Voii/n) - n°\ 

+ P{z n (i,0) < qi + n 1 - f3 V (i/n) - n Q ) 

einf+f 3 <i<bn l '-l 3 

To the first sum in (7.13) we apply a large deviation argument as in (5.29)-(5.30). 
Pick £i, £2 > small enough and take n large enough so that 

(q - ^-^ll^olloo)" 1 {qi + n^Vo^/n) - n e ) < i - e 2 n e 
for all 1 < i < e^e+P. Then 

• i-1 



p{ J> ~ n-"K||oo)*; < qi + n^V^i/n) - n4 
< P^^Xj < i -£ 2 n^| < exp{-m(l -£ 2 n^ -1 )} 



•i=o 

< exp(-Cn 2e r 1 ) < exp(-C 1 n^- fB ) . 

To the last sum in (7.13) we apply Lemma 5.1. Pick < e < g/(u — (3) — 1/2 so 
that e E (0, 1/2). Then for this range of z's 

P{z n (i, 0) < qi + n 1_/3 Vb(i/w) - n Q ) 

< P{z n (i, 0)<qi + n^Voii/n) - ei 1/2+£ } 

< C 2 exp(-C 3 i 2e ) < C 2 exp(-C 3 n 2e(e+/3) ) . 
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Combining the estimates gives 

P{z n (i,0) Kqi + n^PVoii/n)-^} < Cm"- 13 exp(-C 2 n 7 ) 

where 7 > is a new exponent that depends on the earlier constants. The same 
bound is valid for the entire first sum in (7.12). 

By Lemma 5.2, the probability in the second sum in (7.12) is at most 

pjrQraa;] + [2n v qt] -i,n u t) < (nx + 2n v qt - if - C 3 n g 

<exp(-CW 3/2)( ^ /3) ). 

The constant C3 G (0, 1) appeared in front of n e to subsume the difference between 
n"tq 2 + nxq - qi + i 2 /(An v t) in (7.12) and (nx + 2n v qt - if /(An v t) in (7.14). 
Combining the estimates, we get 

J2 P i Z n <r n } <oo. 

n 

Borel-Cantelli now gives (7.11) and completes the proof of Lemma 7.1. □ 



Proof of Theorem, Case 2. Assuming v = 1 + /?, the goal is now 
(7.15) lim n p ~ 1 \z n ([nx\ + [2n v qt],n u t) -n v tq 2 - nxq \ = V(x,t) . 

n^oo I I! 

Let y be a number that achieves the infimum in (2.13) so that V(x,t) = V (y) + 
(x — y) 2 /At. Set i = [ny] in the expression inside the braces on the right-hand side 
of (5.14). Repeat the calculation in (7.2), to get an upper bound for large n: 



(7.16) 



z n ([nx] + [2n u qtln v t) 

< z n ([ny], 0) + r n ([ra/], [nx] - [ny] + [2n v qt\,n v t) 

< nyq + n 1_/ V (y) + n 1/2+£ 

+ n -tq 2 + + nxq - nyq + n^/ 3+£ 

< n v tq 2 + nxq + n 1 ~ f3 V(x, t) + o(n 1_ 



The above steps are justified by Lemmas 5.1 and 5.2, as was done in (7.3)-(7.4) 
below. Again e > needs to be small enough. The estimate n 1//2+e + n( 1+ ^/ 3+e = 
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o(n 1- ^) follows from (3 < 1/2, which itself is a consequence of the assumptions 
v = 1 + (3 and v > 3(3. This gives one half of the goal (7.15). 

For the other half of the proof we can also follow the argument of Case 1. Since 
(3 < 1/2, we can choose g so that 



l + (3 



< 1/2 < q < 1 -(3. 



By the variational formula (2.13), for all i, 

n v tq 2 + nxq + n 1 ~ p V{x 1 1) - 2n e 



< n v tq 2 + nxq + n 



i-P Voii/n)+n i-p ( x 



qi 



n) — n 1 



+ 



4n"t 



— (nx + 2n u qt — i) 



2n g 

2 



Now the argument of Lemma 7.1 can be repeated to conclude that almost surely, 
for large enough n, 

z n ([nx] + [2n v qt],n v t) > n v tq 2 + nxq + 7i 1_/3 "l/(x, t) - 2n Q , 

which together with (7.16) implies (7.15). Case 2 is proved, and thereby Theorem 
4. 
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